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Abstract 

Wc introduce a class of one-dimensional discrete space-discrete time stochastic 
growth models described by a height function ht (x) with corner initialization. We 
prove, with one exception, that the limiting distribution function of h t {x) (suitably 
centered and normalized) equals a Fredholm determinant previously encountered in 
random matrix theory. In particular, in the universal regime of large x and large t the 
limiting distribution is the Fredholm determinant with Airy kernel. In the exceptional 
case, called the critical regime, the limiting distribution seems not to have previously 
occurred. The proofs use the dual RSK algorithm, Gessel's theorem, the Borodin- 
Okounkov identity and a novel, rigorous saddle point analysis. In the fixed x, large 
t regime, we find a Brownian motion representation. This model is equilvalent to the 
Seppalainen- Johansson model. Hence some of our results are not new, but the proofs 
are. 

Key Words: Growth processes, shape fluctuations, limit theorems, digital boiling, random 
matrix theory, Airy kernel, Painleve II, saddle point analysis, invariance principle. 
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1 Introduction 



Growth processes have been extensively studied by mathematicians and physicists for many 
years (see, e.g., 
hansson |2l 
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and references therein), but it was only recently that K. Jo- 
proved that the fluctuations of the limiting shape in a class of growth models 
are described by certain distribution functions first appearing in random matrix theory 
(RMT) 



iOL 41]. Further work by Johansson [29], Prahofer and Spohn [35, gg] and Baik 
and Rains Q strongly suggests the universal nature of these RMT distribution functions. 
These developments are part of the recent activity relating Robinson-Schensted-Knuth 
(RSK) type problems of combinatorial probability to the distribution functions of RMT, 
see e.g. § 1, | |, @, H H H, 0, g| |||, || Q. 



In this paper we analyze a class of one-dimensional discrete space-discrete time stochastic 
growth models, called oriented digital boiling p0| , 21]. Digital boiling dynamics is a cellular 
automaton that models an excitable medium in the presence of persistent random sponta- 
neous excitation. Alternatively, digital boiling models represent contour (constant height) 
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Figure 1: ODB Simulations: The contour lines of ht(x) are drawn in space-time (x,t). In 
these simulations < t, x < 800 and p = 1/2. 



lines for one of the simplest models for growing a connected interface. It is this latter point 
of view we adopt here; that is, we introduce a height function ht{x) that characterizes the 
state of the system. Fig. |l] illustrates the height fluctuations in oriented digital boiling. 
We shall derive various limit theorems for ht(x). We find four limiting regimes: 

1. GUE Universal Regime: x — ► oo, t — ► oo such that p c : = 1 — x/t is fixed and p < p c . 

2. Critical Regime: x — ► oo, t — > oo such that p c : = 1 — x/t is fixed and p ~ p c . 

3. Deterministic Regime: x — > oo, t — > oo such that p c := 1 — x/t is fixed and p > p c . 

4. Finite x GUE Regime: Fixed x and £ — > oo. 

The limit theorems are stated at the beginning of §3. Here is an outline of how they 
are obtained. First we show that ht(x) satisfies a last passage property, i.e. it equals 
the maximum over a certain class of paths in space-time. Then applying the dual RSK 



algorithm [30, |39(], we obtain a reformulation of the problem in terms of Young tableaux. 
This is followed by an application of a theorem of Gessel []l8| (see also p4| ) which gives a 
Toeplitz determinant representation for the distribution function for ht{x). An identity of 



Borodin and Okounkov [1C] expresses the Toeplitz determinant in terms of the Fredholm 
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determinant of an infinite matrix. Finally we use a saddle point analysis (steepest descent) 
to determinine the limiting behavior of the entries, and therefore the Fredholm determinant, 
of the infinite matrix. 

Along the way we identify^] ODB with a first-passage percolation model of Seppalainen [38| 
whose limit law in the universal regime was determined by Johansson [^9|. Thus we could 
have used the analysis in pjj to establish our limit law in the universal regime, or alter- 



natively used Riemann-Hilbert methods J|, 15, 16], to investigate the Toeplitz determinant 



asymptotics. But the method we present is in our opinion more straightforward and tech- 
nically simpler than these, and it is very general. (The Fredholm determinant is easier 
to handle than the Toeplitz determinant, even though they are essentially equal.) Also, 
our analysis permits a nice conceptual understanding of the various limiting regimes. For 
example, the universal regime is characterized by the coalescence of two saddle points; and 
the emergence of the Airy kernel is related to the well-known appearance of Airy functions 
in such a saddle point analysis p^| . 

Even in this simpler approach there are technical details to work out after the saddle 
point analysis gives us the answer. For example in the universal regime we need uniform 
estimates on the entries of the infinite matrix in order to show that the matrix scales in 
trace norm to the Airy kernel. These details are given completely only for this regime. 

In Regime 4 we give an independent proof that the suitably centered and normalized 
ht(x) has a limiting distribution. The proof proceeds through the introduction of a certain 
Brownian motion functional. This leads to some apparently new identities for n-dimensional 
Brownian motion; see ( 4.28| ) below. 

The initial conditions are corner initialization. Due to the fact there is no known sym- 
metry theorem for the dual RSK algorithm |3C, 39 1, we are unable to prove limit theorems 
with different initial conditions, e.g. growth from a flat substrate. From work of Baik and 
Rains |7| and Prahofer and Spohn (3^, |36|, it is natural to conjecture that the limiting dis- 
tribution is now of GOE symmetry and hence given by the analogous distribution function 
in the GOE case Q. 

The table of contents provides a detailed description of the organization of this paper. 



2 Growth Models and Increasing Paths 

In this section we introduce three classes of discrete space and discrete time stochastic 
growth models. Each of these models will have an equivalent path description, but only 
for one of these models are we able to prove limit theorems. Nevertheless, we believe it is 
useful to place this "solvable" case in a larger context. 

We assume that the occupied set of our growth models can be described by a height 
function hf : Z + — » Z + U {— oo, oo}, where Z + is the set of nonnegative integers. Here, time 
t = 0, 1, 2, . . . proceeds in discrete steps. The occupied set at time t is thus given by 

m = {(x, y) e Z + x Z+ : y < h t (x)} . 
A referee points out that this identification can be made at the very beginning. 
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In the models below we use the following one-dimensional neighborhood: (x + jV) = {x — 
l,x} (the oriented case) and assume corner initialization, 

h (x) = { °' iiX = °' (2 1) 

I — oo, otherwise. 

2.1 Oriented Digital Boiling 

The first class of growth rules we call oriented digital boiling (ODB) |^0|, |2l[].Q The rules 
for ODB are 

(1) ht(x) < ht+i(x) for all x and t. 

(2) If h t (x - 1) > h t (x), then h t+ i(x) = h t (x - 1). 

(3) Otherwise, then independently of the other sites and other times, ht+\(x) = ht(x) + 1 
with probability p. (With probability 1 — p, we have ht+\(x) = ht(x).) 

It follows from these rules that for every x and t, ht(x — 1) < ht(x) + 1. 

This process can be readily visualized by imagining the growth proceeding by the addi- 
tion of unit squares starting with the initial square centered at (1/2, —1/2). We denote this 
initial time by placing a in this box. At time t = 1 a box is added to the right (centered 
at (3/2, —1/2)) and with probability p a box is added to the top of the initial box (centered 
at (1/2,1/2)). We place l's in the boxes added at time t = 1. The boxes that are added 
stochastically (Rule (3)) are shaded. An example of this process run for seven time steps 
is shown in Fig. ||. 

2.1.1 Path Description 

As has been observed many times before (see, e.g., [||, [l3|, gl], ^3], [3^, [3S|1 ), a most productive 
way to think about height processes is to introduce the (discrete) backwards lightcone of a 
point (x,t). Precisely, if S = Z + x Z + denotes space-time, then 

Cb(x, t) = {(x , t') G S : < x < x, x < t' < x + 1 — x}. 

For those space-time points in Cb at which a box was added stochastically (according to 
Rule (3)), we place a x. We call such space-time points marked. We define the length of a 
sequence it = {(xi,t±), . . . , (xk,tk)} of (distinct) space-time points in Cb to be k. Such a 
sequence tt is increasing if < xi — Xi-\ < U — ti-% — 1 for i = 2, . . . , k. Let L(x,t) equal the 
length of the longest increasing sequence of marked space-time points in Cb(x, t). (If x < t 
and Cs{x,t) contains no increasing path, then L(x,t) := 0.) For the example in Fig. [2], 
the discrete backwards lightcone £b(3, 7) and an increasing path are shown in Fig. ||[ One 
observes that hj(3) = L(3, 7) = 4. Indeed, this is a general fact. However, before proceeding 
with its proof, it is useful to change slightly the point of view of the process defined by Rules 

2 For spatial dimensions greater than one, visual features of this dynamics resemble bubble formation, 
growth and annihilation in a boiling liquid, hence the process is called digital boiling and oriented refers to 
the choice of neighborhood TV, see Fig. 2 in @ or Feb. 12, 1996 Recipe of 
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Figure 2: Oriented Digital Boiling Process. The number in a box is the time this box was 
added and if the box is colored, then the box was added stochastically according to Rule 3. 
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Figure 3: The backwards lightcone of the point (x, t) = (3, 7) for the process shown in 
Fig. ^. The x's denote the marked points and polyogonal line gives a longest increasing 
path. The length of this path is equal to the number of x 's in the path. This length equals 
h t (x). 
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(l)-(3). Let LT = be a random subset of S to which every point of S belongs with 
probability p. We mark the points of S that belong to II. (Accordingly, we call these points 
marked.) L(x, t) remains the same; namely, the length of the longest increasing sequence of 
marked space-time points in Cb(x, t). With regard to the process, we may intuitively think 
that all the "coins" used in Rule (3) are thrown in advance — of course, many of these are 
ignored as x at time t may become occupied deterministically by Rule (2). Precisely, Rule 
(3) is replaced with 

(3) Otherwise, 

' h t {x) + 1, if (x,t) e n, 



h t+ i(x) 



h t (x), if (a;, t) 011. 



We are now ready to prove the last passage property 
Proposition [21]. h t (x) = L{x,t). 



Proof. We first show that our process is attractive^ in the following sense: Let II and 
II' be two sets of marked points such that LT C LT'. Let ht evolve using LT and h' t using 
LT', then ht < h' t for all t. For if this were not true, then, for some t, h s < h' s , s < t, 
and ht+\{x) > h' t+1 (x) for some x. This, of course, implies that ht(x) = h' t (x). But then 
ht+i(x) = h t (x) + 1 either because of Rule (2); in which case, h' t (x — 1) > h t (x — 1) > 
ht(x) = h' t (x), so (by Rule (2)) h' t+1 (x) = h' t (x) + 1; or, because (x,t) £ LT C II', so again 
h' t+1 (x) = h' t (x) + 1. This is a contradiction. Thus we've established the attractiveness of 
our process. 

The property of attractiveness immediately implies ht(x) > L(x,t), since any increasing 
path of length k will, without the addition of other marked points, cause ht(x) > k. 

We now show that ht(x) < L(x, t). We will show, by induction on k and t, that ht(x) = k 
implies there exists an increasing sequence of marked points of length k in CB(x,t). This 
is obviously true for either t = or k = 0. (Note that ht(x) > means that x < t.) Now 
assume the claim has been demonstrated for all k! < k and t! < t. We can clearly assume 
that h t -\{x) = k — 1, or else we can use the induction hypothesis right away. Therefore, 
we have two possibilities. 

Case 1. ht(x) = ht-\(x) + 1 by application of Rule (2). This means (by Rule (2)) that 
h-t~i( x — 1) = k. Thus by the induction hypothesis, there is an increasing sequence of length 
k in Cb{x — 1, t — 1) C Cb(x, t). 

Case 2. ht(x) = h t _i(x) + l by application of Rule (3'). This means that h t _i(x) = k — 1 
and (x,t — 1) £ II. By the induction hypothesis, Cs(x,t — 1) contains an increasing path 
of length k — 1. Adjoin the marked point (x,t — 1) to the sequence. Observe that the 
increasing property is preserved. This completes the proof of the proposition. 

We summarize this section by noting that ht satisfies for all t > 1, x > 0, 

h t (x) = max{/i t _i(x - l),h t -i{x) + e Xjt } 



where e x ,t = 1 if (x,t) £ LT and otherwise. The initial conditions are (2A). (We take 
ht{— 1) = — oo.) Formulated this way ODB is a "stochastic dynamic programming" prob- 
lem. 



3 For examples of the kind of exotic shapes that can occur from cellular automaton rules without this 
monotonicity property, see WM. 
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2.1.2 The (0, 1)-Matrix Description of ODB 



Without changing the increasing path property, the backwards lightcone Lb of any space- 
time point (x, t) can be deformed into a rectangle of size (t—x) x (x+1). Thus the equivalent 
problem is to fix x and t and to set m = t — x, n = x + 1, and to consider a (0, l)-matrix 
A of size m x n. We number the rows of A starting at the bottom of A and the columns 
of A starting at the left of A. A increasing path in Cb becomes a sequence of l's in A at, 
say, positions ■ ■ ■ , (ikijk)} such that the %i (£ = 1, . . . ,k) are increasing and the ji 

(£ = 1, . . . , k) are weakly increasing. Any such (0, l)-matrix A of size m x n corresponds 
(bijectively) to a twodine array 



w A 



Ji 32 

h 12 



jk 
h 



(2.2) 



where j± < ji < ■ ■ ■ < jk and if jg = jg+i, then ig < ig + \ and the pair (•?) appears in wa if 
and only if the (i,j) entry of A is 1. Note that the upper numbers belong to {1, 2, . . . , n} 
and the lower numbers to {1, 2, . . . , m}. For example, the matrix 

/000 1 00 1\ 

11110 11 

1 1 1 

10 110 11 

1 1 

V 1 1 / 

maps to the twodine array 



.4 



w A 



11122233344456667777 
34524523535613451356 



(Recall the convention for row labels.) As an example, a longest increasing path (of length 
5) is indicated in bold typeface. We remark that one can compute the length of an increasing 
path by patience sorting 0] on the bottom row of wa (from left to right) with the rule that 
a number is placed on the left most pile such that it is less than or equal to the number 
showing in the pile. Patience sorting on the above example results in the five piles 



5 
5 

6 6. 



If N denotes the number of l's in a random m x n (0, l)-matrix A, then the above 
mappings imply that for any nonnegative integer h, 



Prob(h t (x) < h) 



Prob (h t (x) <h\N = k) Prob (N = k) 

fc>0 



4 We have chosen both a nonstandard labeling of A and a nonstandard bijection A 
increasing path property remains (essentially) the same under the bijections. 



wa so that our 
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£ )/(l -p) mn - fe Prob(L m>nife < h) (2.3) 



fc=0 



where L mn ^ is the length of the longest increasing path in a random (0,l)-matrix A with 
l's (or equivalently, in the associated wa)- 



2.1.3 Tableaux Description of ODB 

The dual RSK algorithm [|30l , |3l| is a bijection between (0, l)-matrices A of size m x n and 
pairs (P, Q) such that P* (the transpose of P) and Q are semistandard Young tableaux 
(SSYTs) with sh(P) = sh(Q) where the elements of P are from {1,2, ... ,m} and the 
elements of Q are from {1,2,... ,n}. In terms of the associated wa, (|2.2| ), one forms P 
by successive row bumping of the second row of wa starting with %\ and with the rule an 
element i bumps the leftmost element > i. Thus each row of P is strictly increasing. A 
fundamental property of the dual RSK algorithm is that the length of the longest strictly 
increasing subsequence of the second row of wa equals the number of boxes in the first row 
of P. 

If d\(M) denotes the number of SSTYs of shape A with entries coming from {1, 2, . . . , M}, 
then the number of pairs (P, Q) of fixed shape A in the above dual RSK algorithm is 

d\>(m)d\(n) 

where A' is the conjugate partition. (Conjugate since P* is a SSYT.) Since there are (™ n ) 
(0, l)-matrices with k l's, 

Prob (£ m , n ,fc <h) = — £ d y (m)d x (n) = — ^ d x (m)d y (n). 

V k / Ahfe \ k I Ahfc 

A x <h l(\)<h 

And hence from ( |2.3| ) 

mn 

Prob(h t (x)<h) = (l-p) mn Y / r k £ d A (m)d A '(n) 

fc=0 Al-k 
<(A)<h 

where r = p/(l — p). Observe that for |A| > mn, d\(m)d\i(n) = 0. (A SSYT with 
entries from {1,2,..., M} can have at most M rows.) If V denotes the set of all partitions 
(including the empty partition), then the above sum can be summed over all partitions 
without changing its value, 

Prob(fct(s) <h) = (l-p) mn r\ x \d x (m)d x ,(n). (2.4) 

xev 
e(x)<h 

Comparing (p.4|) with Johansson's Krawtchouck ensemble results establishes the equiv- 
alence of ODB with the Seppalainen-Johansson model. 
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2.1.4 Application of Gessel's Theorem and the Borodin-Okounkov Identity 

Gessel's theorem |l8|, 44] is 

r^s x (x)s x (y) = D h (<f) 



l(\)<h 



where s\ are the Schur functions (see, e.g. |39tl ) and Dh(<p) is the hxh Toeplitz determinant^ 
with symbol 

oo 00 

<p(z) = IIU - x 3 z)~ l - Vjrz- 1 )- 1 . 

3=1 j=l 

If we apply to both sides of this identity the automorphism u (see Stanley Iffil , pg. 332), 
u>(s\) = s\i, to the symmetric functions in the x- variables we obtain 



J2 r\ x \ Sx/ (x)s x (y)=D h ( i p) (2.5) 

l{\)<h 

where now the symbol is 

oo oo 

ip{z) = \{{l + x jZ ) H(l- yj rz- l )-\ (2.6) 

3=1 3=1 

Recalling the specialization ps* (see Stanley p9|], pg. 303), we apply ps^ to the x- variables 
and ps^ to the y-variables in Gessel's identity ( |2.5| ) and observe^] that the resulting LHS is 
precisely the RHS of (2.4). Since the specialization ps^ is a ring homomorphism, we may 



apply it directly to the symbol (|2.q) . Doing so we obtain 

Prob (h t (x) <h) = (l- p) mn D h (<p) (2.7) 

where 

(p(z) = (l + z) n {l-r/z)~ m . (2.8) 
This derivation required r < 1. However, by (|2.3| ) the left side is a rational function of r, 



and analytic continuation shows that (2.7) holds for all r > if in the integral representing 
the Fourier coefficients of <p the contour has r on the inside. 

The Borodin-Okounkov identity expresses a Toeplitz determinant in terms of a Fred- 
holm determinant of an infinite matrix which in turn is a product of two Hankel matrices. 
Subsequent simplifications of the proof by Basor and Widom || extended the identity to 
block Toeplitz determinants. We now apply this identity to the Toeplitz determinant fl2.7p . 
First we find the Wiener-Hopf factorization of <p(z): 



5 If is a function on the unit circle with Fourier coefficients (f>k then T n (4>) denotes the Toeplitz matrix 
6t-3')»,j=o,...,»-i and D„((f>) its determinant. 
6 Note that ps^ s\ = d\{n) which follows from the combinatorial definition of the Schur function. 
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where 

<p+{z) = (l + zr, ^(z) = (l-r/z)- m . 
Define K h acting on £ 2 ({0, 1, • • •}) by 

oo 

K h (j,k) = J2( l P-/<f+)h+j+e+i (<p+/<p-)- h -k-l-i- (2-9) 
1=0 

The Borodin-Okounkov identity is then 

D h (<p) = Z det (I- K h ). 

Since the determinant on the right tends to 1 as h — > oo as does Prob (ht(x) < h), we have 
Z = (1 — p)~ mn . Thus we have derived a representation of the distribution function of the 
random variable ht(x) in terms of a Fredholm determinant, 

Prob (h t (x) <h)= det (I - K h ) . (2.10) 

This derivation also required r < 1. As above, analytic continuation shows that (|2.10|) holds 
for all r > if in the integral representing the Fourier coefficients of y?_ /ip+ the contour has 
r on the inside and —1 on the outside. (In fact the contour must have —1 on the outside 
no matter what r is.) 

A somewhat different direction (and one we do not follow here) is to apply isomonodromy 
and Riemann-Hilbert methods |], [ll], |27j] directly to the Toeplitz determinant D^ip). This 
would result in the identification of D^ip) as a T-function of an integrable ODE. 



2.2 Inhomogeneous ODB 



In ODB the probability p appearing in Rule (3) is independent of the site x. Inhomogeneous 
ODB replaces Rule (3), for each site x € Z+, with 

(3 X ) Otherwise, then independently of the other sites and other times, h t +i(x) = h t (x) + l 
with probability < p x < 1 and ht+i(x) = h t {x) with probability q x := 1 — p x . 

Since the dual RSK algorithm is a bijection between (0, l)-matrices A and pairs (P, Q) 
such that P t and Q are SSYTs with col(^4) = type(P) and row(^4) = type(Q) |3^], we have 



VToh{h t (x)<h)=af---q™ dx(m)s x ,(r) 



(2.11) 



l{X)<h 



where, as before, m = t — x, but now r = (ro, . . . , r x , 0, . . .) with rj := Pj/qj- The proof of 
( 2.11] ) is straightforward and similar to the proof of the analogous result in [26]; therefore, 
we omit it. The right hand side of ( |2.11 ) clearly reduces to (2.4) in the homogeneous case. 
We again apply Gessel's theorem to obtain the Toeplitz determinant representation 



Prob (ht(x) < h) 



q? 



whereQ 

<p(z) = {l-l/z)- m f[(l + 

3=0 

7 The homogeneous case of (2.12) does not directly reduce to (|2.8| 



It does after z 



(2.12) 

z/r which corre- 



sponds to a similarity transformation of the Toeplitz matrix. 
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Application of the Borodin-Okounkov identity results in a Fredholm determinant represen- 
tation for this distribution function. Observe that from either (|2.11| ) or ( |2.12| ) it follows 
that Prob(h t (x) < h) is a symmetric function of {po,pi, ■ ■ ■ ,p x )- This property opens the 
possibility for an analysis of the spin glass version of ODB which we plan to address in 
future work. 

2.3 Weak ODB and Strict ODB 

Here are two natural variants of the ODB. We let the "spontaneous increase" in Rule (3) 
apply after Rule (2) has already taken effect to get weak ODB: 

(1') ht(x) < ht + \{x) for all space-time points (x,t). 

(2') If h t (x — 1) > ht(x), then ht(x) = h t (x — 1) else h t (x) = ht-i(x). (Here h t is an 
intermediate height function.) 

(3') Independently of the other sites and other times, ht+i(x) = ht(x) + 1 with probability 
p. (With probability 1—p, ht+i(x) = ht(x).) 

In strict ODB we require that the left neighbor is restec^ for the spontaneous increase. 
(We take ht(x) = — oo for x < which in this model implies ht(0) = for every t.) 

(1") ht(x) < ht+i(x) for all space-time points (x,t). 

(2") If ht(x - 1) > h t (x), then h t+1 (x) = h t (x - 1). 

(3") Otherwise, if a? — 1 is rested at time t, ht{x — 1) = ht(x) then independently of 
other sites and times, ht+x(x) = ht(x) + 1 with probability p (ht+i(x) = ht(x) with 
probability 1—p.) 

In a similar way one shows 

• In weak ODB, ht(x) equals, in distribution, the longest sequence of positions 
in a random (0, l)-matrix of size m x n (m = t — x + 1, n = x + \) which have entry 
1 such that in are ji are both weakly increasing. (The lower left corner of the matrix 
is fixed to be a 0.) 

• In strict ODB, ht(x) equals, in distribution, the longest sequence (ii,je) of positions 
in a random (0, l)-matrix of size m x n (m = t — x, n = x) which have entry 1 such 
that ig are jg are both strictly increasing. 

3 Limit Theorems 

In this section we derive limit theorems for the distribution function Prob(/tt(x) < h) for 
ODB. Our starting point will be the Fredholm determinant representation ( ^.10|) . This 
distribution function is a function of four variables, x, t, h and p; and accordingly, there 
are several asymptotic regimes: 

8 The height at a site cannot increase at two consecutive times, i.e. it must rest for one time unit before 
it is allowed to increase. 
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Figure 4: In the GUE Universal Regime, the left figure displays the limiting shape, c\, as 
a function of x/t and the right figure displays the normalization constant, c 2 , as a function 
of x/t. In both cases, p = 1/2. 



(1) GUE Universal Regime: Let x — > oo, t — > oo such that p c := 1 — x/t < 1 is fixed. For 
fixed p < p c define 



ci := 2p c p - p + 2\f pp c (l -p)(l -p c ), 
c 2 : = (p c (l-p c )) 1/6 (p(l-p)) 1 / 2 



1 + 



'(l-p)(l-Pc) 



We will show that 



Prob 



h t (x) - cit 
c 2 t 1 / 3 



PPc 



< S 



Pc 



I- Pc V 1 - P 



v 



(3.1) 

2/3 

'(3.2) 



F2(s) 



where ||40 | 



^2(5) = det — ^Airy) = ex P J ( x ~ s )l( x ) 2 d% 



(3.3) 



Here ^Airy * s ^ ne °P era tor with Airy kernel acting on L ((s,oo)) (see ( |3.5[) below) 
and g is the (unique) solution of the Painleve II equation 

q" = sq + 2q 3 

with boundary condition q(s) ~ Ai(s) as s — ► 00. The limiting shape, c±, and the 
normalization constant, c 2 , as functions of x/t are shown in Fig. |I] for p = 1/2. The 
probability density, f 2 = dF 2 /ds, is shown in Fig. ||. 

(2) Critical Regime: Let x — > 00, t — > 00 such that 

x = (l-p)t + o(VS). 
For fixed A £ Z + we will show that 

Prob (h t (x) - (t - x) < -A) 
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. 5 r 




Figure 5: The density /2(a) = dF2/ds where F2 is defined by ( |3.3| ). The distribution 
function F2 has mean \i = —1.77109, standard deviation a = 0.9018, skewness S = 0.2241 
and excess kurtosis K = 0.0935. 



converges to a A x A determinant. One can think of this as 

(3) Deterministic Regime: For x — > 00, t — ► 00 and fixed p > p c , we will show that 

Prob (ht(x) = p c t) — ► 1. 

(4) Finite x GUE Regime: Fix x and let i — ► 00, then we will show that 

V(p(l-p)0 1/2 / 

converges to the distribution of the largest eigenvalue in the GUE of (x + 1) x {x + 1) 
hermitian matrices, denoted below by F^+i E ■ 



3.1 GUE Universal Regime 

It is convenient to use the variables m = t — x and n = x + 1 rather than x and t and to 
translate back to the space-time variables at the end. We assume p < p c := m/(n + m). 
(This is asymptotically 1 — x/t as defined above.) Further, when there is no chance of 
confusion, we denote the random variable hf(x) by H. (We reserve lower case h to denote 
the values of H.) Set h = cm + sm 1 / 3 , where c will be determined shortly, and a = n/m. 



(In this notation the condition p < p c is ar < 1.) For any f the matrix ((- 



K h (j,k)) 



has the same Fredholm determinant (the determinant of / minus the matrix) as (Kh(j, k)). 
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We shall show that for a particular v and a certain constant b > this matrix scales to a 
kernel with the same Fredholm determinant as 



^Airy(sM36) 1/3 + x, s/v^b) 1 ' 3 + y), (3.4) 

on (0, oo), where 

;>oo 

K Mry (s + x, s + y) = Ai(t + s + x) Ai(t + s + y) dt. (3.5) 
J o 

This gives 

lim Prob ( H ~, C ™ <s]= F 2 (s/v(3b) 1/3 ). (3.6) 

Here is what we mean by scaling. Any matrix (M(j, k)) acting on £ 2 (Z + ) has the same 
Fredholm determinant as the kernel M([x], [y]) on L 2 (0, oo) and this in turn has the same 
Fredholm determinant as M m (x, y) = m 1 ' 3 MQm 1 / 3 x], [to 1 / 3 ?/]). If this kernel has the 
limit k(x, y) we say that the matrix (M(j,k)) has, after the scaling j — > m 1 / 3 x, k — > 
m 1 / 3 y, the limit k(x, y). If M m (x, y) converges to k(x, y) in trace norm then the Fred- 
holm determinant of (M(j,k)) converges to that of k(x, y). And if (M(j,k)) were the 
product of two matrices each having scaling limits in Hilbert-Schmidt norm (under the 
same scaling, of course), then the Fredholm determinant of the product converges to the 
Fredholm determinant of the product of the limits. This is what we shall show in our case. 

There is a slightly awkward notational problem. Since h is always an integer and 
h = cm + sm 1//3 , the quantity s as it appears here and the analysis which follows is 
not completely arbitrary. What we actually show is that if h and m tend to infinity, and s 
is defined in terms of them by the formula h = cm + sm 1 / 3 , then 

Prob (H < h) — F 2 (s/w(36) 1/3 ) (3.7) 

uniformly for s lying in a bounded set. From this we easily deduce ( |3.6| ) for fixed s, which 
now has a different meaning. These observations are important when one tries to estimate 
errors. It can be shown that the difference in ( |3.7| ) is 0(?n~ 2 / 3 ). But the difference between 
the right side of (|3.6j ) and the probability on the left can only be expected to be 0(m^ 1 ^). 
The reason is that if the quantity s' is defined by cm + s'm 1 / 3 = [cm + sm 1 / 3 ] then the 
probability is within 0(m _2/3 ) of F 2 (s' /v(3b) 1/3 ), but s-s' is very likely of the order m ' . 



3.1.1 The Saddle Point Method 

The matrix (Kh(j,k)) is the product of two matrices, the matrix on the right having j, k 
entry (c^ + /92-)-h-j-fc-i an d the one on the left having j, k entry [ip^ / ip + )h + j + k + i. Notice 
that the first vanishes if h + j + k + 1 > m so we may assume that all our indices j and k 
satisfy h + j + k < m. We have 

(^ + /^)^_,_ fe _ 1 = ±J(l + z) n (z-r) m z- m+h+i+k dz (3.8) 
= (-l) fe +J+ fc J- J(\ + z ) n (r - z) m (- z )~ m + h +i+ k dz, 
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and a similar formula holds for If we set 

tj){z) = (1 + z) n (r - z) m (-z)- {1 - c)m 

then 

(-l)fc+i+*( ¥ , + /y,_)_^ i _ Jfe _ 1 = 2, | ^(z)(-z)- 1/3 +^ fc dz 

and 

(-i)^-+ fc ^_/^ + ), +i+fc+1 = JLy ^r 1 (-«)-™ i/8 -i-*- a (fe. 

The contours for the first integral surrounds while the contour for the second integral 
has r on the inside and —1 on the outside. The restriction h + j + k < m is the same as 
sm i/3 + j + < (l _ c ) m . if we make the replacements j — ► m 1 / 3 x, — ► m 1 / 3 y these 
become 

— f ij){z) (-z) ml/3 ( s+x+ ^ dz, — f ^(z)- 1 (-z)~ ml/3 (-+*+iH-2"»- 1/3 ) dz. 

Our restrictions become s + x + y<(l — c) m, 2 / 3 . For convenience we replace s + x + y by 
x, and we want to evaluate 

^ | ^(z) (-z) ml/3 *dz, ^{z)- 1 (-z)- ml/3 *~ a (3.9) 

asymptotically. Our restriction is now x < (1 — c) m 2 / 3 . 
To do a steepest descent we have to find the zeros of 

d , , . . n m (1 — c) m 

— log ^ (*) = T - — + " — , 

dz \ + z z — r z 

or equivalently the zeros of 

(c + a) z 2 + (c + r — cr — a r) z + r (1 — c) . 

(Recall that a = n/m.) The discriminant of this quadratic equals zero when 

1 



1 + r 



(2y/ar +(l-a)r). (3.10) 



This is the value of c we take|] The critical probability is the condition c = 1, i.e. p c 
m/(m + n). The single zero of the quadratic is then at u = —v where 

(1 — r) c + (1 — a) r 1 — \Jar 
2(c + q) ~ 1 + Vo^' 



If there were two critical points, or if we took the negative square root in (3.10), the Fredholm deter- 
minant would tend exponentially to either zero or one. It is only for this value of c that we get a nontrivial 
limit. 



16 



u r 

Figure 6: The steepest descent curves as described in the text. 



Note that < v < 1 since < p < p c . (It is because u < that we used powers of — z 
rather than z in the definition of ip.) We write 



1 d 3 

6b := — — t log ip(z) 
m az a 



2a 



+ 



2(1 - c) 



(1 + u) 3 (u — r) 3 
2(v^ + y/Tf 



r y/a (1 + r) 3 (1 — y/ar) 
The quantity 6 is positive since ar < 1. In the neighborhood of 2 = u, 

V(z) ~^(-u)e mb( ^ )3 . 



(3.11) 



The steepest descent curves will come into u at angles ±7r/3 and ±27r/3. Call the former 
C + and the latter C~ . For the integral involving ip{z) we want 1^(2)1 to have a maximum 
at that point of the curve and for the integral involving ^(z) -1 we want (V'(-z)l to have a 
minimim there. Since b > the curve for i^{z) must be C + and the curve for ^(z)" 1 must 
be C~ . Both contours will be described downward near u. The curve C + will loop around 
the origin and close at r, the upper and lower parts making an angle there depending on 
c while C~ will loop around on both sides and go to infinity with slopes depending on c. 
(That have these forms follows from the fact that the contours cannot cross and, since 
the only critical point is at z = u, the contours can end only where ip, respectively ip" 1 , is 
zero.) The steepest descent curves are shown in Fig. ^. 

Proceeding formally now, consider the ip{z) integral and make the substitution z — > 
u + z = —v + z. Then the old — z becomes the new v (1 - z/v) ~ve~ z / v , and recall (|3~Tl| ). 
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If we make these replacements in the integral we get 

1 /•ooe~ i7r / 3 
M u ) v ™ 1/3 x / ^f-mWxz/v dz 

2-Kl Jooe in / 3 

The contour can be deformed to the imaginary axis since we only pass through regions 
where Kz 3 is negative. If we then set z = —iC,/m 1 / 3 the above becomes 

2n 7_oo (So) 1 / 6 
If we recall that x was a replacement for s + x + y we see that the matrix with j, k entry 

has the scaling limit 

1 Ai((s + x + y)/v(3b) 1 / 3 ). 



Similarly the matrix with j, k entry 

(-l) h +^(u)v ml/3s+ 3 +k (^/tp + ) h+3+k+1 

has 1/v 2 times exactly same scaling limit. Hence the scaling limit of u fc_J Kh(j,k)), which 
has the same Fredholm determinant as (K^j, k)), is the product of these scaling limits, 



1 



v 2 (36) 2 /3 J 



Ai((s + x + t)/v{3b) 1/3 ) Ai((s + 1 + y)/v(3b) 1/3 ) dt 



1 — - K Airy ((s + x)/v(3b) 1 / 3 , (s + y)/v{3b) 1 ' 3 ). 



v{3bY 

And, as promised,^] this kernel has the same Fredholm determinant as 

K Ahy {s/v{3b) 1 ' 3 + x, s/v(3b)^ 3 + y). 

3.1.2 Convergence Proof 

Now for the justification. We have to obtain not only the pointwise limit, but uniform esti- 
mates to establish convergence of the operators in trace norm. We first obtain asymptotics 
under the assumption that x lies in a bounded set. (Notice that x > s always.) We begin 
with 

t/;(z)z ml/3x dz, (3.12) 



2-7TZ Jci 

and denote by the portion of C + which lies within e of the critical point u c >. 

Lemma 1. If in ( p.!2[ ) we integrate only over the error incurred is 0(\if)(u) \ e~ Sm ) for 
some 5 > 0. 



10 T he t ime c onst ant ci = p c c and the normalization constant C2 = p X J Z v (36) 1//3 . A computation then 
gives Q and Q . 
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Proof. Define 

a(z) = — log ij)(z) = a log (1 + z) + log (r — z) + (c — 1) log (— z). 
m 

Its maximum on C + (it is real- valued there) is cr(u) and it is stricly less than this on 
the complement of C e in C + . Therefore i(j(z) /i/}(u) = 0(e~ Sm ) for some 5 > on the 
complement while z ml 3x = e°( ml 3 ). This gives the statement of the lemma. 

Lemma 2. We have as m — > oo 

m 1 / 3 MuY 1 v ~ ml/3x — [ M z ) z ml/3 x dz = -(36)- 1/3 Ai(x/w(36) 1/3 ) + Ofm" 1 / 3 ) 

2m J 

uniformly for bounded x. 
Proof. Near z = u 

a(z) = a{u) + b{z - uf + 0{{z - uf)), 
log(— z) = \ogv + — (z — u) + 0((z — u) 2 ) = \ogv — -(z — u) + 0((z — u) 2 ). 

U V 



Hence, using Lemma 1, we have 



— / ip(z) z ml/3 x dz 
2m Jc+ 



= 0(\l(>(u)\e~ Sm ) + 1p(u) v ml/3x -}— f e mb(z-u) 3 -m 1 ^x(z-u)/v+0(m(z~u) 4 +m 1 ^ 3 x(z-u) 2 ) ^ 

2m Jc+ 

We show that removing the O term in the exponential in the integrand leads to an error 
0(m~ 2 / 3 ) in the integral. This error equals 

mb(z— u) 3 — m 1 / 3 x(z— u)/v f ^p{m{z—u) A +m 1 / 3 x{z—u) 2 ) _ -y \ ^ 
em b(z-u) 3 - m ^ 3 x(z-u)/v+0( m (z-ur+m^ 3 x(z-uf) Q ,, _ ,4 + _ 

ct 

Now the exponential has the form 

rnb(l+rn)(z—u) 3 —m 1 / 3 x(l+ri2)(z—u) /v 
& j 

where the rji can be made arbitrarily small by taking e small enough. If we make the 
substitution z — u = ^/m 1 / 3 the error becomes 

2/3 [ f> Kl+ m )C 3 -x(l+ m )C/v 



The integral is now taken over a long contour lying in thin angles around the rays |arg£| = 
7r/3, with ends having absolute value at least a constant times m. This integral is clearly 
bounded, uniformly in m for bounded x. 

Therefore with the stated error we may remove the O terms from the exponential in 
the original integral. Then we make the same substitution. The integrand is exponentially 
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small at the ends of the resulting contour. Therefore if we complete it so that it goes to 
infinity in the two directions ±7r/3 the error incurred will be exponentially small. 
We have shown that 

m 1 / 3 Mu)~ l w -™ 1/3 * _L / M z ) z ml/3 x dz = / e K3 - x < /v d( + 0(m" 1/3 ). 

2iri J 2m Jooe™/ 3 

If we deform the contour to the imaginary axis and make the substitution £ — ► — iQ then 
the last integral, with its factor, becomes 

_ 1 r e ib^+ixc/v d( - = _(36)-i/3Ai( a ./ u (36)V3) > 
2tt J-oo 

This proves the lemma. 

The second integral in ( |3.9| ) is similar. 

Lemma 3. We have as m — > oo 

i 1 / 3 V(n) v ml/3x — [ ^{z)- 1 z- ml/3 x dz = -(36)- 1/3 Ai(xM36) 1/3 ) + CKm" 1 / 3 ) 



m 

2iri jq- 
uniformly for bounded x. 

Proof. The derivation is essentially the same. The exponentials are replaced by their 
negatives and the directions ±7r/3 are replaced by ±2tt/3. The fact that C~ is unbounded 
causes no difficulty since the integrand now behaves at infinity like a large negative power 
of z. We get the same Airy function in the end, as we have already seen. 

Now for the tricky part. We need estimates that are uniform for all x and where the 
error term contains a factor which is very small for large x. In fact we shall show that the 
statements of Lemmas 2 and 3 hold, uniformly for all x, when the error terms are replaced 
x in the exponential can be improved to a constant times x 3 / 2 but that 
makes no difference.) To do this we have to be more careful and use the steepest descent 
curves for the full integrands in fl3,9| ), not just for the factors tf) . We consider in detail 
only the first integral in (|3.9[) ; as before, the second is treated analogously. 

Set 

if)(z, c') = (1 + z) n (r - z) m (- z )-( x - c ')™ = i/;( z ) 
We are interested in the asymptotics of 

I(c') = ^-Ji;(z,c')dz (3.13) 

when d — c = to -2 / 3 x. Our condition on x says that d < 1 and, in view of what we already 
know, we may assume x is positive and bounded away from zero, so d > c. 

We let C be the steepest descent curve for ip(z,d). This curve now passes vertically 
through one of the critical points of ip(z, d). For d > c there are two critical points 



, -(1 - r)d - (l-a)r±\/ ((1 + r) d + (a - 1) rf -4ar 
uf, = 



2{a + d) 
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which are real and satisfy 



1 < u c , < -v < uj, < 0. 



To determine which critical point our curve passes through we consider the function 

a(z, c) = — log ip(z, d) = a(z) + (c - c) log (—z). 
m 

The critical points are the zeros of a z (u c >, d). (Subscripts here and below denote deriva- 
tives in the usual way.) We use the fact that uj; are smooth functions of 7 = \/c' — c and 
compute, recalling that cr z (u, c) = cr zz (u, c) = and observing that dd jdrf = when 7 = 0, 



d 



du. 



± 



■j-cr zz (u^,c) = a zzz (u,c)^- . (3.14) 

7=0 



dj 7=0 «7 

The first factor on the right is positive (we denoted it by 66), while 

du^, 



= ±(3 (3.15) 

7=0 



where 







(ar) 1 /^! + r )3/2 



(v^ + v^) 2 ' 

Since a zz (u^, d) = when 7 = we deduce that for d close to, but greater than, c we have 

Vzz(u+,,c) > 0, cr zz (u~,,c) < 0. 

These inequalities hold for all d since the second derivative can be zero only when d = c. 
This shows that the steepest descent curve C for tp(z, d) passes through u^,, because on the 
curve \t/j(z,d)\ has a maximum at the critical point. (Similarly the steepest descent curve 
for ijj{z, d)~ l passes through u~,.) To make the notation less awkward we write u c > instead 
of ut. First, we have the analogues of Lemmas 2 and 3. 

Lemma 4. Given e > there exists a 5 > such that 1(d) = 0(\ij)(u, d)\e~ Sm ) if d — c > e. 

Proof. The function a(z,d) is decreasing for u < z < u c > since it decreases near and 
to the left of u c > and has no critical point in this interval. Hence a(u c >,d) < a(u,d) 
so cr(u c ',d) — a(u,d) is negative and bounded away from zero for d > c + e. Since the 
maximum of \ip(z, d)\ on C is at z = u c i the statement follows. 

In view of Lemma 4 we may assume in what follows that d — c is as small as we please. 
We denote by C e the portion of C which lies within e of the critical point u c i. 



Lemma 5. If in the integral (3.13), in which we integrate over C, we integrate only over 
C e the error incurred is 0(\ip(u, d)\ e~ Sm ) for some 5 > 0. 

Proof. The maximum of \i)(z,d)\ l / m onC occurs at u c > and it is stricly smaller on the com- 
plement of C £ in C than it is at u c < . Therefore the integral in question is 0(\ip(u c i, d) \ e~ Sm ) 
for some 5 > 0. Since a(u c >,d) < a(u,d), as we saw in the proof of the last lemma, this 
one is established. 
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Because of Lemmas 4 and 5 we need only compute the behavior of ip(z, d), or equivalently 
a(z, d), for z near u c i . Recall that u = u c = —v. 

Lemma 6. We have 

(i)a(u c ,,d) = a(u) + (d-c)logv-^(d-cf/ 2 + o((d-c) 2 ) ; 

(ii) a z (u c >,c') = ; 
(Hi) a zz (u c r,c) = 66/3 \fd — c + 0(d — c) ; 
(iv) a zzz (u c >,d) = 6b + 0(Vd - c) . 



Proof. From ( |3.15 ) and the fact that u c i is a smooth function of 7 (or directly) we see that 

u d - u = f3y/d - c + 0(c - c). (3.16) 



Consequently, since u = —v, 

^L = 1 -H^r^ + +0 (c' - c ). (3.17) 

U V 

Now since a z (u c /, d) = we have 

d d u 

—a(u c >,c) = —a(z,c) ^ ^ = log(-iv) = log u + log . 

Integrating with respect to d from c to d and usi ng (|3.17| ) we obtain (i). Of course (ii) is 
immediate. As for (hi) and (iv), these follow from fl3.14| ) and ( 3.15| ) and the fact that u c > is 
a smooth functions of 7. 

Lemma 7. The conclusions of Lemmas 2 and 3 hold uniformly for all x when the error 
terms are replaced by 0(e~ Sm ) + 0(m~ 1 ^ 3 e~ x ) for some 5 > 0. 

Proof. We consider ( |3.12|) , which is 1(d) with d — c = m~ 2 l 3 x. Putting together Lemmas 
5 and 6 we deduce that 

1(d) = 0(\ip(u) v {c '~ c)m \ e~ Sm ) + ip(u) v {c '- c)m e -^ {c '- c)3/2m x 

1 f e rnb{z-u c ,)' i +^mb(5V(^{z-u c ,)' 2 +0{rn[(c'-c) 2 + {c'-c)\z-u cl \ 2 + ^ 

If d — c = m~ 2 l 3 x the exponential factor equals e~~^ x3/ while the integral equals 

1 f e mb(z-u c ,) s +3m 2 ^' i bx 1 ^ 2 p(z-u c ,) 2 +0(m^ 1 /' A x 2 +m 1 ^ 3 x\z-u c ,\ 2 +m 2 ^ 3 x 1 / 2 \z-u c ,\ 3 +rn\z-u c ,\ 4 ]) ^ 

2-ki Jc e 

Now C e , rather than looking like two rays near the critical point, looks like one branch of 
a hyperbola. 

Note that by Lemma 4 we may assume that d — c = m~ 2 l 3 x is as small as desired. It 
follows that the exponent, without the 0(m~ l / i x 2 ) term, can be written 

m(b + m) (z - u d f + 3m 2/3 (b + r/ 2 )x 1/2 /? (z - u c >) 2 , 
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where, if e is chosen small enough, the r\i can be made as small as desired. Upon making 
the variable change z — u c > = f^/m 1 / 3 the integral becomes 



m 1/3 / e (b+vi)( 3 +3(b+v2 )x^pe d( - 
2-ki 



taken over a long contour in the right half-plane on which |arg£| > tt/3 — 77, with another 
small rj. The integral here is uniformly bounded. 

To take care of the term 0(m _1//3 x 2 ) in the exponential in the original integral, observe 
that if m~ 2 / 3 x is small enough then m~ l ^x 2 will be at most a small constant times x 3 / 2 
and so 

Thus removing the term from the exponential leads to an eventual error 

That removing the other O terms from the exponential leads to the same error is seen as 

it was in the proof of Lemma 2 — the substitution in the integral representating the error 

results in an extra factor m -1 / 3 and there is the exponential factor e~^ x3/ outside the 

integral. 

After removing all the O terms and making the variable change z — u c > = Q/m 1 /' 3 the 
integral becomes 

m " V3 f e K 3 +3bx 1 / 2 K 2 



2vri 

taken over a long contour in the right half-plane on which |arg £| > tt/3 — rj. Completing the 
contour so that it goes to infinity in the directions arg£ = ±7r/3 leads to an exponentially 
small error. It follows that (the first part of) the lemma holds with the negative of the Airy 
function in the statement replaced by 

If we complete the cube and make the substitution £ — ► £ — fix 1 / 2 this becomes, upon noting 
that 3b (3 2 = l/v, 

ooe - " r / 3 

JP-xC/vq = _(36)-V3 Ai (a;/(„(36)V3)). 

ooe l7r / 3 

The second part of the lemma is analogous, just as the proof of Lemma 3 was analogous 
to the proof of Lemma 2. 

We have now shown that if we set j = m 1 / 3 x, k = m lj/3 y then 

(_l)^+;+fc m l/3 ^ 1/3 ^+^( 93+ / (/9 _)_ ?i _ J _ fc _ 1 _> -(36)- 1 /3 Ai((s +x + y)/ V ^ h )V^ 

and the difference between the two is 0(m~ 1 / 3 e~( x+y ^) + 0(e _<5m ). It follows easily from 
this that if we denote the matrix on the left, without the factor m 1 / 3 , by (M(j,k)) and 
the kernel on the right by A(x,y) then the kernel m 1//3 M([m l l 3 x], [m 1 / 3 y]) converges in 
Hilbert-Schmidt norm to the kernel A(x,y) on (0, oo). (Recall that j and k are at most 
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0(m). Therefore the error term 0(e~ ) can only contribute an exponentially small error 
to the norm and so can be ignored. Similarly we can let our indices j and k run to infinity.) 
Thus, under the scaling j — > m 1//3 x, k — > m 1//3 y the matrices with j, k entry 

scale in Hilbert-Schmidt norm to the kernel A(x, y). Similarly so do the matrices with j, k 
entry 

Therefore the product of the matrices scale in trace norm to the (operator) square of the 
kernel, which is the Airy kernel ( |3.4| ). This or completeness the justification. 



3.2 Critical Regime: p ~ p c 

When p = p c (pir = the analysis of the previous section must be modified. We set 
h = m — Ah (Ah = 0, 1, 2, . . .) and introduce the new ip 

^ = (l + z ) n ( z - r ) m 

and the corresponding new a 

a(z) = — logt/> = alog(l + z) + \og(z — r). 
m 

The saddle point now occurs at z = with c"(0) = — a(l + a). Thus in the neighborhood 
of z = 

^(z) ~ (-l) m r m e - ma{1+a) z2/2 . (3.18) 

Since ((p+ / (f-)-h-k-j-i vanishes for h+j+k+1 > m, we can again assume h+j + k < m 
which becomes the condition j + k < Ah. As before our starting point is the integral 
expression 

( V >+/<p-)- h . j . k - 1 = ^-J i>(z)z- m+h+ j +k dz 

where the contour is a circle centered at with radius p < 1. Taking this p sufficiently small 
so that we may use the approximation ( 3.18 ) on the integrand, we obtain after making the 
change of variables 

1/2 

c 



{T^L) ' z = z /s, 



( (_ 1 ^m r m sj+ k-Ah+i t^L if Ah - j - k - 1 = 2L = 0, 2, 4 . . . 

if Ah — j — k — 1 = odd integer. 

11 See the remark at end of this section. 



24 



Our second integral is 

(<p-/<p + ) h+j+k+1 = ^-J ^{zT X z m - h -i- k - 2 dz 

where the contour has —1 on the outside and r on the inside. We deform the contour to 
the imaginary axis going from ioo to — ioo with an infinitesimal indentation going around 
to the left. The part of the contour lying in the right half plane is exponentially small 
because of the factor (1 + z)~ n and can therefore be neglected. For the integral along the 
imaginary axis we can replace ip by Q3.18 ) with an error that is exponentially small. Thus 
the above integral is asymptotically equal to 

_— m /•—too „ „ 

{-l)™ r — / e z2 / s2 z Ah ^ k ~ 2 dz, 

which in turn equals 

1 P — OO „ 

/ -Qm t -Afc-.7-fc-l r -m gAfe-j-fc-l 1 / e ~( 2 ^Ah-j-k-2 ^ 

2iri Joo 

where there is an indentation above £ = 0. If we now substitute £ = y/i, the above integral 
becomes 

1 r0 + 

t 1 \m^A/t-,7-fc-l r -m gAh-j-k-l 1 / g-i < (A/i-j-fc-l)/2-l ^ 

4vri ioo 

The contour starts at +00, loops around in the positive direction and then returns to 
+00. This last integral is Hankel's integral representation of the T function. Thus 



We now use these two asymptotic expressions along with the condition j + k < Ah in 
|2.9Q to obtain (after a short calculation) 



K h m~ { -^P- ± \^l^-j)T{i + k -^-y (3.19) 



When £ + (k — j)/2 is a nonpositive integer, the product of the sine and gamma functions 
is replaced by 

The factor {—S) k ~^ may be dropped when computing the determinant det(/ — K^) since it 
does not change its value. We evaluate this determinant and display the results for Ah < 9 
in Table |. 

Remark. Since m and n are integers it is extremely unlikely that p = p c = m/(m + n). 
If p is irrational this never occurs. However the preceding analysis shows that if ar = 
1 + o(m -1 / 2 ) rather than 1 then in the integrals one gets extra factors (1 + z ) ( ml ^ 2 \ 
Then after the substitution z = S( this drops out since S = Oim" 1 / 2 ). The upshot 
is that the asymptotics hold for any p when m and n go to infinity in such a way that 
m/(m + n) = p + o(m -1 / 2 ). 
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Ah 


lim„ woo Prob (H — m < — Ah) 


Numerical Value 





1 




1 







1 


l 

2 







■5 




2 


1 
4 


1 

2tt 


9 


AO A £ 1 

U8451 x 


1 n — 2 
10 


3 


1 

8 


3 

8tt 


5 


63379 x 


10~ 3 


4 


1 

16 


i_ 1 29 
' 37T 1 96tt 


1 


17616 x 


10" 4 


5 


1 

32 


41 145 
r 144tt 2 768tt 


8 


22908 x 


io- 7 


6 


1 

64 " 


32 | 1169 1249 
1357r :i 1 3840?r 2 102407T 


1 


92570 x 


lO" 9 


7 


1 

128 


49 | 198827 8743 
225^ 1 9216007T 2 122880tt 


1 


50565 x 


io- 12 


8 


1 

256 


4096 10289 , 5773487 145603 
' 23625ir 4 360007T 3 1 34406400tt 2 3440640tt 


3 


92048 x 


10 -16 


9 


1 

512 


15376 5528469 , 279234531 436809 
~ r 91875T 4 250880007r a 1 2569011200ir 2 18350080tt 


3 


42524 x 


1Q -20 



Table 1: Limiting Distribution when p ~ p c 



3.3 Deterministic Regime: p > p c 
3.3.1 Large Deviations Approach 

Assume that p > p c . Then there exists an e > so that n/m approaches (1 + e)(l/p — 1). 
To simplify the statements, we will just assume that n = (1 + e)(l/p — l)m. 

Imagine the random m x n matrix A from §2.1 as the lower left corner of an infinite 
matrix of 0's and l's, created by the independent coin flips. Fix a position > 1) 

in this infinite random matrix. Define J as the column index of the first entry, from left to 
right, with a 1 on the row above and in the columns larger or equal j. Then define 
(iu j) = J — j- In the example given, £(3,1) = and £(5,1) = 3. 

Now create a sequence of i.i.d. random variables £i,£2 5 ---j as follows. Let £1 equal 
the column index minus one of the first 1 on the first row. Then let £2 = £(i.i+£i)> £3 = 
£(2,i+£i+£ 2 )> • • •• T ne basic observation is that, since we are always taking the best positioned 
1 on the next line, we have equality of the two events 

{there is an increasing path of length m in A} = {£1 + . . . + £ m < n}. 

Therefore, we need to show that 

Prob(£i + . . . + U > n) 

goes to exponentially as m — > 00. However, Prob(^i = i) = p{l — p)" 1 , i = 0, 1, . . . and so 
= 1/p — 1. By elementary large deviations (e.g. §1.9 in ]T7|), 

-m" 1 log P(£t + . . . + U > n) -> 7(e) 

where an elementary calculation shows 

7(e) = (1/p - 1)(1 + e)log(l + e) -p _1 (l + e - ep) log(l + e - ep) 

which is positive whenever e > 0. 
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3.3.2 Saddle Point Approach 



For completeness, we show how the saddle point method gives the same result. Thus we 
show that when p > p c (or ar > 1) 

det(J - Kh) — > 

exponentially as m —> oo even when h = m — 1, thus establishing assertion 1(c) in §3 with 
exponential approach to the limit. 

As we saw at the beginning in the last section we need only consider the entries Kh(j, k) 
when h+j + k < m, which in the present situation means j = k = 0. Our claim is therefore 
that K m -i(0, 0) — ► 1 exponentially asm-> oo. The first integral to consider is 

(<p+/<P-)- m = ^~J iiz) z~ l dz = m = (-r) m . 
The second integral is 

(¥>-/V+)m = — ■ J lj){z)~ 1 Z~ 1 dz. 

Recall that the contour here surrounds and has —1 on the outside, r on the inside. The 
critical point for steepest descent is at z = u where 

a 1 ar — 1 

+ = 0, 



1 + u u — r a + 1 

The steepest descent curve will pass vertically through this point and go to oo in two 
directions. But notice that since u is positive, in order to deform our original contour to 
this one we have to pass through z = 0. The residue of the integrand there equals (— r)~ m 
and so 



{<p-l<p+) m = (-r)- m + —J Mz)-'z-' dz, 

where now the integral is taken over the steepest descent curve. This integral is asymp- 
totically a constant times m -1 / 2 times the value of the integrand at z = u, and this value 
equals (— l) m times 

' a(r + l) \- am / r + 1 
a + 1 / \a + \, 

Our claim is therefore equivalent to the statement that this is exponentially smaller than 
r~ m , which in turn is equivalent to the inequality 

n a 

( r + 1)°+! _ > r 

[ + ' (a + l)«+i 

It is an elementary exercise that this is true for all r > except for r = 1/a, when equality 
holds. But in our case r > 1/a so the inequality holds. 
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3.4 Finite GUE Regime: Fixed x and t — > oo 
3.4.1 Saddle Point Calculation 

We return to fl3.8|) and this time set 

h = m + sm l l 2 = pm + s m 1 / 2 , 

1 + r 

and make the substitutions j — ► xm 1 / 2 , fc — > ym}' 2 to write the integral (13.81) as 



^~J( 1 + Z T ( r ~ Z T (-z)- m/{1+r) (-z)( s+x+y ^ ml/2 dz. (3.20) 

Now we set 

^(z) = (r- z) m (-z)" m /( 1+r ), 

which is the main part of the integrand. There is a single critical point, z = —1, and at 
this point d 2 /dz 2 \ogip{z) is equal to 

r 

m Tr l — \2 = ™>p(l-p). 
(1 + r) z 

This is positive and so the steepest descent curve is vertical at the critical point; it goes 
around the origin and closes at z = r. The main contribution to the integral comes from 
the immediate neighborhood of the critical point. If we make the variable change 

.= -1 + C 



'771 

and take into account the other factors in the integrand we see the integral is asymptotically 

_ ( r + 1 ) m r (J_) n (3. 21) 
2iri i^ioo \\/m Jm 



Now 

1 r io ° r, p -b 2 /4a 

2m J-ioo 2Jam 



and so 



2iri J-^ 2y/cm db n ^r(2 v / a) n+1 n \2sfa J ' 

(H n are the Hermite polynomials.) Hence our first integral ( |3.20| ) is asymptotically equal 
to — (r + l) m /y / m n+1 times this expression with 

a = -^p(l-p), b = s + x + y. 
Thus we have shown that the matrix with j, k entry 

(_l)-*-i-*( r + i)-™ m »/2 ( ¥?+ / ¥J _)_ /l _ i _ fc _ 1 
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scales to the operator on (0, oo) with kernel 

1 (s+x+y) 2 /4a f s + x + y 

v^(2^) n+1 H 2^E 

with a as given above. 

Next, with the same substitutions in the integral, 

(j?-/ V+ ) h+j+M = -L [(l + z)->(z-r)- m z m ->'->- l -' 2 dz 



2ni 



The contour here encloses and r and has —1 on the outside. The steepest descent curve 
should go through the critical point — 1 horizontally. We deform the given contour to a curve 
starting at — oo + Oi, going above the the real axis, looping around z = — 1 clockwise, then 
back below the real axis to — oo — Oi. The original contour can be deformed to this because 
the integrand is small at oo. The main contribution is again in the neighborhood of z = — 1. 
Making the same variable change as before leads to an integral which is asymptotically 

_ ( r + 1 )~ m f ( _C_V n e -^ P (i-p)( 2 +(s+x+yK (3.22) 
2iri J \\/mJ \fm 



where now the contour is a circle going around £ = counterclockwise. Using now the fact 



1 r „ 1)/2 „ Jn— 1 

— f C~ n e~< +b< d£ = - — e~ b ' ia - e -(C-6/2v^) 2 

2m J s s (n-1)! dC"" 1 

Q (n-l)/2 / ft 

H n -i 



C=o 



(n-1)! " _1 V2V5, 
we find that the matrix with j, k entry 

(_l)"+i+*( r + l) m rn n ' 2 (<p_/<p + ) h+j+k+1 

scales to the operator on (0, oo) with kernel 

aC"- 1 )/ 2 f s + x + y 



(n-1)! 2^ 

Combining, we see that the product of the two matrices (aside from a factor (— 
which does not affect the determinant) has scaling limit the operator with kernel 



s + z + y 



1 [°° e -(s +X+ zf/Aa H ( S + X + Z \ 

^2^a(n-l)\J n \ 2^E ) n ~ X \ 2^E 



dz. 



Instead of a direct evaluation of this last integral, we will not evaluate our £ integrals 
( 3.21| ) and ( 3.22 ), but rather consider them as integrals with variables Ci an d C2, combine 
and integrate with respect to z. We see that the scaled kernel for the product is 



i Jo" 1 1 (r) " e a ^-®- {s+x+z)(l+is+z+y)(2 dz dCi d( 2 , 
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where the Ci contour is a vertical line described upward and the C2 contour goes around 
counterclockwise. If the vertical line is to the right of the circle we can integrate first with 
respect to z, yielding 



1 



2 



^ 2 J J vc 2 ; C1-C2' 

Let's call this L n (x, y). This is when n = 0, and 

L k (x, y) - L k ^{x, y) = ~f [ ^^<MM^)Ci+<*+tO& d(l d ^ 

s2 



4vr2 J J C: 



This integral is a product and we can use the computations we did above to see that it 
equals 

1 -(*+*)* /4a rr ( s + x \ ( s + y \ 

2*^WE{k - 1)! \ 2^E J k ~ x \2^Ej- 

If ipk are the oscillator wave functions^ then this equals 



times the factor 



1 (s + x\ fs + y\ 



e -{s+xf/8a e ( s+y )2/8a_ 



It follows that if Kjj.n is the Hermite kernel then 

L n (x, y) = e~^ 2 ^ * K H , n S -±E-) e^^. (3.23) 

We deduce that 

lim Prob (H < pm + s m 1 ^ 2 ) 

m— >oo 

is equal to the Fredholm determinant of 



1 „ f s + x s + y \ 



over (0, 00), or equivalently the Fredholm determinant of Kh,ti(x, y) over (s/2-^/a, 00). It 
is notationally convenient to introduce 

a 2 := 2a = p(l — p) 

and to define 

F^ UE (s) := lim Prob ( < s) . 



This equals the Fredholm determinant of Kh,u over ( S /V2, 00) and is equal to the distri- 
bution of the largest eigenvalue in the finite n GUE.p| 

12 The oscillator wave functions are <fik{x) '■= e~ x 1/2 Hk(x) /y/2 k k\ -k 1 / 2 and form an orthonormal basis for 
L 2 ((0,oo)). The Hermite kernel is KH,n(x,y) := J22=o <Pk(x)<fk(y)- 

13 Our normalization of F^ UE differs from the usual one |33, 142] by a factor of y/2, i.e. the usual normal- 
ization is the Fredholm determinant of the Hermite kernel over (s, 00). 
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3.4.2 Moments of F® UE 



From the theory of random matrices, e.g. 42], we know that the distribution function 
det(J — Kh,ti) has an alternative representation as an n x n determinant. Explicitly, 



? GUE 



(s) = det 6ij 



/V2 



(pi(x)ifj(x) dx 



0<i,j<n-l 



where ipj are the oscillator functions previously introduced. This last representation implies 
that the F^ UE are expressible in terms of elementary functions and the error function with 
increasing complexity for increasing values of n. In the simplest case, n = 1, F^ UE is the 
standard normal; a result easily anticipated from the original formulation of the growth 
model. The next simplest case is n = 2, 



zpGUE 



\ 1 1 " 

(s) = 4"^ e 



1 



s e 



2 3 /2 0f 

The moments of F„ E are, of course 



* 2 /2 + I (1 _^ se -, 



2vr 



erf(s/V2) + - erf^/v^) 2 . 



/oo 
S^f^ UE (s)ds,j = 1,2,... 
-oo 



where f£ UE = dF^ UE /ds. For 1 < n < 5 we have, 
First Moments: 



Mi(l) = 
Ml (2) = 

Mi (3) = 

Mi (4) = 

Mi (5) = 

Second Moments: 

M2(l) 
M2(2) 

M2(3) 
M2(4) 

M2(5) 



0. 



-= « 1.128379, 

a/VT 



27 
8^ 



1.904140, 



7 



+ 



475 



+ 



475 



48^2 vr 3 / 2 1280F 64 vr 3 / 2 

13715 16975 41145 

+ 



arcsin(l/3) » 2.528113, 



4O960r 41472^2 vr 3 / 2 2048 vr 3 / 2 



1. 

2, 

9\/3 

3 + -f- » 4.240490, 



arcsin(l/3) ps 3.063268. 



4 + 



47T 

16 



6.940420, 



r 155\/5 
5 - — „ + 



2495 



+ 



499 



864 7T 2 108^/3 vr 54v/3vr 2 



arcsin(l/4) rj 1.977575. 
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n 


A 4 


Approx 


a 1 


Approx 


S 


K 


2 


1.12838 


1.251 


0.72676 


0.645 


0.08465 


0.01053 


3 


1.90414 


1.989 


0.61474 


0.564 


0.11862 


0.02192 


4 


2.52811 


2.594 


0.54907 


0.512 


0.13749 


0.03042 


5 


3.06327 


3.118 


0.50426 


0.476 


0.14972 


0.03683 


6 


3.53861 


3.585 


0.47101 


0.448 


0.15838 


0.04184 


7 


3.97026 


4.011 


0.44497 


0.425 


0.16490 


0.04586 


8 


4.36822 


4.405 


0.42379 


0.407 


0.17001 


0.04917 


9 


4.73920 


4.772 


0.40609 


0.391 


0.17414 


0.05195 



Table 2: The mean (//) and the variance (a 2 ) of H™, 2 < n < 9, {H™ has distribution 
function F^ UE ) are compared with the approximations ( [3.24 ) and ( |3.25 ), respectively. Also 
displayed are the skewness (S) and excess kurtosis (K) of . H°° has S ~ 0.2241 and 
K ^ 0.0935. 



Third Moments: 



Ma(l) = 0, 

u 3 (2) = 4= ~ 3.949327, 

\/7T 



297 

/i 3 (3) = — ^= « 10.472769, 

16V7T 

333 7109 7109 

u 3 (4) = = — + -7TT + ^ arcsin(l/3) ss 20.378309, 

32^2 7T 3 / 2 256ttV2 128 vr 3 / 2 w 7 

2595475 259385 778155 

M3(5 = = — - + = + ^777 fn 33.432221. 

82944^^/2 8192^F 4096 ^ 3 / 2 



Fourth Moments: 

/i 4 (l) = 3, 
/i 4 (2) = 9, 

/i 4 (3) = 19 + » 28.096927, 

2 7T 

496 

/x 4 (4) = 33 + — =- » 63.384348, 
3V37T 

, 4 (5) = 51 + ™*4 + -gg°P- + " 5 I 5 , arcsin(l/4) «, 117.872208. 
P4V 7 1296^ 2 324^ 162^^2 V / ' 

Let denote the weak limit m — > oo, n fixed, of 

H — pm 

and the weak limit m — > oo, n — ► oo, a = n/m fixed, of 

1 



u(36)V3 m i/3 
32 



(ff — cm) . 



Figure 7: The densities fn UE > 1 < n < 7. For increasing values of n the maximum of fn UE 
moves to the right. 

(Thus the distribution functions of iT^° and H°° are F^ UE and F%, respectively.) For 
a — > 0, c = p + 2oyfa + O(a), (36) 1 / 3 ~ er/a 1 / 6 , and t> ~ 1. Proceeding heuristically, 

H ~ cm + m 1/3 v(3&) 1/3 #°° 

~ pm + 2cr v / a"i + m 1/3 cra" 1/6 if 00 

~ jim + am 1 / 2 ^2y/n + —jygj • 



Thus we expect 



and hence 



/6 



£(tf£°) « 2v^+^p, E(H°°) = —1.77109 . . . , (3.24) 
Var(iJ^) « Var(g°°) y ar /g-oo\ _ q.8132 (3.25) 

For 2 < n < 9, these approximations are compared with the exact moments in Table [2|. We 
also compute the skewness and the excess kurtosisp^ of F^ UE . 
The densities f^ UE , 1 < n < 7, are graphed in Fig. ^. 

4 Brownian Motion Representation in the Finite x GUE 
Regime 

Let B(t) = (Bo(t), . . . , B x {t)) be the (x + l)-dimensional Brownian motion. Let F be the 
following functional on continuous functions / = (fo,---,f x ) from [0,1] to R 21 , which 

14 The skewness of a random variable X is E (( X ~ M ) 3 ) and the excess kurtosis is i? (( A "~ M ) 4 ) — 3. Here 
fj, = and cr 2 = Var(X). 
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satisfy /(O) = 0, 

F(f) := max{/ (i ) + fi(h) - fi(t ) + ... + f x {t x ) - f x {t x ^) : < t < h < . . . < t x = 1} . 
Note that F is continuous in the L°° metric. Finally, let 

M x := F(B). 

Theorem. For x 6 Z+, and t — > oo, we have 

h t (x)-pt d 



-^M x , 



where a 2 = p(l — p). 

Proof. First, the path representation tells us that we change ht{x) by at most a constant 
if we only obey the increasing property within the same column. That is, we change L(x, t) 
to L'(x,t), where L'(x,t) is the longest path (xi,ti),i = l,...,k, of marked points such 
that < ti — ti-i — 1 if Xi = Xi-i, while < Xi — rrj_i < ti — — 1 if Xi 7^ x^±. Thus, 
the first observation is 

\L(x,t) — L'(x,t)\ < x. 

Let SI equal the length of the longest increasing sequence of points (t, i), < t < k. 
Then 

L'(x, t) = max{^ + - + . . . + £%, - S^_ t : < h < k 2 < . . . < k x = t - x}. (4.26) 

Now for every fixed i, S\ is independent of S 3 k for j / i. Let X l k equal the indicator of the 
event that (i, k) is a marked point. Of course, S l k = J2e=i X\. 

Let S 1 (t), t G R_|_, equal S l k when r = k and be obtained by linear interpolation off 
the integers. Moreover, let S l be the centered versions S l {r) = S 1 (t) — pr, and S(t) = 
(5 (r), . . . , S x (t)). For < r < 1, define 

Xt{T) - = 77r' 

then the standard invariance principle (see, e.g. p7| , Ch. 7.) implies that X t converges as 
t — > 00 in distribution to the (x + l)-dimensional Brownian motion B. 
Now define 

L"(x,t) = max{5°(to)+5 1 (ti)-5 1 (t )+. • .+S x (t x )-S x (t x ^) : < t < h < . . . < t x = t}, 

(4.27) 

then 

\L'(x,t) -pt- L"(x,t)\ < 5x. 

(The linear interpolation gives an error of at most four at each ti and we incur an additional 
x by replacing t — x by t.) Note now that (by making a substitution t[ = tj/t) 



L"(x,t) 

^7r- F(Xi) - 
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It follows (by continuity of F) that the theorem holds with L"(x,t) in place of h t (x), but 
this is clearly enough. 

Remark 1. We should note that this theorem clearly holds in more general circumstances. 
For example, we could make every x count an independent random number of jumps. We 
would get the same theorem, with the only assumption that the said random number has 
finite variance. 



Remark 2. The theorem also holds for random words over an alphabet with n letters [ 44 1 , 
except that the x + 1 Brownian motions are not independent, but they have to sum to 0, 
so the covariances equals (n — l)/n 2 when i = j and —1/n 2 otherwise. In the case of 
two equiprobable letters, the limiting distribution of the centered and normalized length of 
the longest weakly increasing subsequence in a random word is equal to the distribution of 
the random variable 

X = m^Boit) + B l (l) - Brit)) 
= 2max(Bo(t))- J Bo(l) 
= 2M - N. 

(M denotes the random variable maxo<t<i Bo(t) and N denotes the random variable Bq(1).) 
From the reflection principle it follows (see, e.g. pg. 395 in [|l^|) that the joint density of 
(M, N) is 

hi N(m, n) = J— (2m - n) e ~( 2m ~ n ) 2 / 2 , for m > 0, m > n. 

V 7T 

Thus the density of X equals^ 



x 10 o 



fx(x) = / fM,N(m, 2m - x) dm = J- x 2 e x 12 . 

JO V 7T 

Remark 3. Limiting distribution of the centered and normalized ht(l): Here we have 
M a = max (iM*) + (B 2 (l) - B 2 (t))) = max (J3i(i) - B 2 (t)) + B 2 (l) = M + N. 

(The random variables M and N are defined by the last equality; and therefore, are not to 
be confused with the random variables of the previous remark.) Note that N is standard 
normal. Since (B\ — B 2 )/V2 is the standard Brownian motion, M equals, in distribution, 
y/2\N\ again by the reflection principle. Even though M and N are not independent, 
E(Mi) = V2E(\N\) = 2/^/tt. Moreover, the conditional distribution of N given the entire 
path of W := B\(t) — B 2 (t), < t < 1, depends only on its final point VF(1). Given this 
final point S equals s, the distribution is normal with mean — s/2 and variance 1/2. That 
is, if J~t is the Brownian filtration for W, S = W(l), then 

Prob (N e dn\Fi) = Fmb{N £ dn\S = s) = - e - [x+s/2)2 dn. 

7T 

lj C. Grinstead, in unpublished notes, also found a random walk interpretation of the two-letter random 
word problem and used this to determine the limiting distribution in this case. 
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This makes it immediately possible to compute the second moment of M\, since 

E(MN) =E(E(MN\F±)) = E(ME(N\F l )) = -E(MS)/2 = -E ({M '/ 'V2)(S '/ '\/2)) = -1/2, 

by a straightforward computation with the joint density above. Therefore E(M 2 ) = 
E{M 2 ) + E{N 2 ) + 2E(MN) = 2. 
In this way, the density of M\ is 

f Ml (x) = E (Prob (M 1 = x\M,S)) 

= J f M , s (m, s) Prob (Mi = x\M = m,S = s) dmds 

roc rm 1 1 

= dm -f M:N (m/V2,n/V2)—e-^- m+n ^dn. 

JO J-oo 2 WTT 



An explicit evaluation shows this last integral equals, as it must, f% (x). 

Remark 4- Since the distribution function of M x equals FS^, it follows from RMT |?3| 
that we have the alternative representation 

Prob (M x <s)=Cn ( ■■■ f A(x) 2 e"^ x ? dx x ■ ■ ■ dx n (4.28) 

J— oo J — oo 

where 

A(x) = A(xx, . . . ,x n ) = Yi ( x i~ x j) 

l<i<j<n 

is the Vandermonde determinant, c" 1 = 1!2! • • • n\ (27r)™/ 2 , and n = x + 1 > 2. In the 
context of Brownian motion, can one directly prove ( [4,28j )? 

Remark 5. For connections between Brownian motion exit times and random matrices, 



see Grabiner [19| 



Remark 6. The Brownian motion functional M x has appeared previously in Glynn and 
Whitt [25], and consequently ( 4.28| ) provides an exact formula for the limiting distribution 



of the departure time of the first (x + 1) customers from n single server queues. Glynn and 
Whitt also consider the case when x = t a , < a < 1, and prove what would, in our setting, 
be the following limit theorem 

h t (x)-pt M x 
lim = — = a := nm — —, 

t^oo aJtx :r - >00 Jx 



with both limits in probability. They conjectured that a = 2, and this was later proved by 
Seppalainen [37] via a hydrodynamic limit for simple exclusion. We note that our paper 



proves that a = 2 as well, by a completely different route. Namely, one only needs to apply 
the result (see, e.g. [Q]) that the largest eigenvalue in the finite n GUE scales as 2-y/n.^ 
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